Departure from the conical dispersion in epitaxial graphene 
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The TT bands of epitaxiahy grown graphene are studied by using high resolution angle resolved 
photoemission spectroscopy. Clear deviations from the conical dispersion expected for massless 
Dirac fermions and an anomalous increase of the scattering rate are observed in the vicinity of the 
Dirac point energy. Possible explanations for such anomalies are discussed in terms of many-body 
interactions and the opening of a gap. We present detailed experimental evidences in support of 
the gap scenario. This finding reveals a fundamental intrinsic property of epitaxial graphene and 
demonstrates the possibility of engineering the band gap in epitaxial graphene. 



With the continuing demand for even smaller transis- 
tors, we will soon reach the physical limits of silicon as 
a host material. This has motivated the search for novel 
materials to replace silicon in the past decade. Among 
the various candidates, carbon-based materials and in 
particular carbon nanotubes have shown impressive po- 
tential for nanoelectronics applications [T . Most recently 
attention has turned to planar graphene sheets as an even 
better alternative thanks to their remarkable electronic 
transport properties at room temperature [2[ [3[ [H [5l E] , 
such as an electron mobility that is even higher than that 
of silicon. Although the possibility of using graphene 
to make transistors has been recently explored both in 
epitaxially grown graphene [5] and in graphene ribbons 
ill ISI |9j , the biggest hurdle to overcome is that graphene 
is a zero gap semiconductor. Over the past few years, 
different ways to induce a gap in the electronic spectra 
and generate a finite mass for the Dirac fermions have 
been proposed. These range from engineering graphene 
nanoribbons [TOl |TT] or quantum dots [T2j |13] , where the 
gap is the result of the quantum confinement of the elec- 
trons in these structures [ij, to applying a voltage to 
bilayer graphene, i.e. by making the top and bottom lay- 
ers inequivalent [T5j [161 Ull EH] , or through the growth 
of graphene on a substrate which will break the carbon 
sublattice symmetry [T9l. 

Here we present a detailed study of the electronic struc- 
ture of single layer epitaxial graphene with the main focus 
on the TT bands by using angle resolved photoemission 
spectroscopy (ARPES). The dispersion of the tt bands 
shows overall agreement with conical dispersions of Dirac 
fermions. However, clear deviations from this conical dis- 
persion are observed near the Dirac point. We examine 
the two possible explanations for the experimental data: 
many-body interactions and the opening of a gap at the 
Dirac point. Our data suggests that the gap opening is 
a more likely scenario. 

The monolayer graphene film was grown epitaxially on 
the (0001) surface of an n-type 6H-SiC by thermal decom- 
position of SiC, as previously reported [18, 20, 21, 22l. 
The surface oxide of SiC is removed either by anneal- 
ing in situ at ~ 900C under silicon flux or by hydrogen 
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FIG. 1: (a) LEED pattern of single layer epitaxial graphene 
on SiC taken at 180eV. (b) Schematic drawing of the structure 
of single layer graphene on top of the buffer layer and the SiC 
substrate, (c) Intensity map as a function of in-plane and 
out-of-plane momenta (/cyand kz) at -1 eV. 



etching before being loaded into the ultra high vacuum 
chamber. The substrate is then heated by electron beam 
to 1400C and monitored by low-energy electron diffrac- 
tion (LEED) until graphene is detected. The substrate 
preparation and annealing conditions are crucial knobs 
for determining the quality of the sample as well as its 
thickness [181 1211 123] . We have measured several samples 
that have been grown under different conditions. Al- 
though the width of ARPES peaks can be affected by 
sample quality, the effects discussed in this paper are 
always present in different samples. The sample thick- 
ness was determined by Auger, x-ray photoemission spec- 
troscopy (XPS) and angle resolved photoemission spec- 
troscopy (ARPES), which consistently showed one mono- 
layer thickness for the samples studied here. ARPES 
measurements were carried out at Beamline 12.0.1 of the 
Advanced Light Source (ALS) in Lawrence Berkeley Na- 
tional Lab with an energy resolution better than 25 meV. 
The samples were measured at 25K. 

Fig. 1(a) shows a typical LEED pattern of the as-grown 
monolayer graphene at 180eV electron energy. The 1x1 
spots for SiC are clearly observed together with the 
(6a/3 X 6y/3)R30 reconstruction pattern due to multiple 
diffraction between the SiC substrate and the monocrys- 
talline graphene overlayer. The (6a/3 x 6\/3)R30 is an 
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FIG. 2: (a-c) ARPES intensity maps taken at Ef, -0.4 eV and 
-1.2 eV respectively on single layer graphene. The dotted line 
shows the Brillouin zone of graphene. (d) schematic drawing 
of the dispersion in single layer graphene and the relative 
energies for data shown in panels a-c. (e) Dispersion of single 
layer graphene measured along a high symmetric direction 
through the K point (see black line in the inset). 



indication of the presence of the carbon rich layer that 
develops before graphene can be grown [24 , see Fig. 1(b). 
This carbon rich layer, also known as the 'buffer layer', 
shows the same a bands as graphene [25 , suggesting that 
these two have the same lattice structure. However, the tt 
bands are absent because the tt orbitals have hybridized 
with the dangling bonds from the substrate [26 . The 
buffer layer is electronically inactive near E^? and thus 
one single layer of graphene plus a buffer layer is neces- 
sary to give rise to the tt band dispersions. In addition to 
Auger and XPS measurements, we have determined the 
number of layers by measuring the band structure with 
ARPES by measuring the dispersions along kz- This is 
a very effective method as the resulting number of quan- 
tum states are directly associated with the number of 
graphene layers [I8l [23l J7 . The absence of tt band 
dispersion along the kz direction in two Brillouin zones 
(BZs) is a strong indication that the sample studied here 
is a single layer graphene. 

Figs. 2 (a-c) show ARPES intensity maps at constant 
energies of E^, -0.4 eV and -1.2 eV. At E^ (Fig.2(a)), 
the intensity map shows a finite and almost circular con- 
tour centered around the K point. As the binding energy 
increases, the contour first decreases in size and becomes 
a point at -0.4 eV (Fig.2(b)). Beyond -0.4 eV, the size of 
the contour expands again. This behavior is consistent 
with the presence of Dirac fermions, where a conical dis- 
persion centered at the K point is expected. The Dirac 
point energy is shifted to 0.4 eV below E^, which shows 
that the as-grown graphene is electron-doped [23, 27]. 
At higher binding energy (Fig.2(c)) the high intensity 
region in the intensity map deviates from the circular 
shape. Similar trigonal distortions have been reported 
for graphite Note that the intensity along the cir- 

cular contour is not isotropic and is strongly suppressed 



on the left (first BZ) for energy above the Dirac point 
energy (Fig. 2 (a)) and on the right (second BZ) for en- 
ergy below the Dirac point energy (Fig. 2(c)). This is a 
well known effect in graphite and is due to the ARPES 
dipole matrix element which suppresses or enhances the 
intensity in different BZs [29] , 

Fig. 2(e) shows the overall dispersion of the tt bands 
measured for a symmetric cut through the K point. Fol- 
lowing the maximum in the intensity plot, two cones dis- 
persing in opposite directions (one upward and another 
downward) can be clearly distinguished, in overall agree- 
ment with the presence of Dirac fermions. The Dirac 
point energy, defined as the midpoint between the va- 
lence and conduction bands at the K point, occurs at 
-0.4 eV since the sample is electron doped, as already 
discussed above. More importantly, in addition to this 
shift of the chemical potential and hence E/^, deviations 
from a conical dispersion are observed in the vicinity of 
the Dirac point. The first, hardly visible on the energy 
scale of this figure, is at approximately -0.2 eV below 
E F and is due to coupling to phonons [30l |311 [32] . The 
second occurs near E/) and is the main focus of this pa- 
per. In particular, we observe that the valence band and 
the conduction band do not merge at a single point at 
El), but instead the top of the conduction bands and the 
bottom of the valence bands stop before E^? and at the 
K point there is finite intensity over an extended energy 
region between these two bands (see shaded area in the 
figure) . 

To gain further insights into the nature of the observed 
anomaly, it is important to perform a quantitative anal- 
ysis of the data. This is usually done by extracting the 
peak positions and the peak width by fitting the mo- 
mentum distribution curves (MDCs) at fixed energy and 
the energy distribution curves (EDCs) at fixed momen- 
tum. By fitting the MDCs using Lorentzian functions, 
the dispersion can be extracted from the peak positions 
and the scattering rate, a quantity proportional to the 
imaginary part of the self energy ImS, can be obtained 
from the peak width. However, one of the caveats of the 
MDC analysis is that it is based on the assumption that 
whenever there is a peak in an MDC there are electronic 
states. 

In Fig. 3, we present data taken along two lines (panel 
a): one through the K point (panel b), and another one 
parallel to it but far from the K point (panel c). These 
particular cuts are convenient because the intensity is 
strongly suppressed on one side of the K point, thus al- 
lowing the MDC to be fitted with single peaks. This 
allows for a more reliable fit of the data. The cut far 
away from the K point is considered because it definitely 
has a gap, due to the conical nature of the dispersion. 
Performing an MDC analysis first on the cut through 
the K point, the anomalous region of Fig. 2 is manifested 
through a kink in the dispersion and a sudden increase 
of the scattering rate (pointed to by arrows in panels c 
and d). Such anomalies might be considered to be due 
to self-energy effects. An appealing explanation is that a 
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FIG. 3: (a) Schematic drawing (not to scale) for the cuts shown in (b) and (c) in the BZ of graphene and the conical dispersion. 
(b,c) Dispersions measured through and off the K point. The dotted white and dark gray lines are dispersions extracted from 
the MDCs. (d) MDC width as a function of energy. 
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FIG. 4: (a,b) Simulation of conical dispersions with a gap of 
150 and 400 meV. (c) Extracted MDC width as a function of 
energy for data shown in panels a and b. 



decay through plasmon emission is responsible for the de- 
viation from the conical dispersion, as recently proposed 
[33| 134) l35] . However, performing a similar analysis on 
the cut far away from the K point reveals similar features. 
In both cases, we can identify a region between the con- 
duction and valence bands where the intensity decreases 
(see horizontal arrows), the dispersion deviates from the 
linear behavior (see dashed line in panel b-c), and the 
scattering rate shows a sudden increase (see panel d). 
These striking similarities cast doubt on the validity of 
the many-body interaction scenario [33], as this should 
be able to account for very similar features far away from 
the K point where a gap is certainly present and is the 
natural explanation of the data. We therefore propose 
that these similarities are simply a manifestation of a 
gap opening at the Dirac point. 

It is therefore misleading to discuss the dispersion and 
scattering rate in the gap region and in general to perform 
an MDC fit near the bottom or top of a band, despite 
the possibility that MDC peaks might be present inside 
the gap region (panel c). Such peaks exist in our data 



due to the finite width of the EDCs at the bottom of the 
conduction band and the top of the valence band. This 
can be better illustrated in the simulated data shown in 
Fig. 4. Dirac-like dispersions with finite band gaps of 150 
meV (Fig. 4(a)) and 400 meV (Fig. 4(b)) and hence finite 
effective masses are used. The MDC width is modeled to 
scale linearly with energy with finite constant broaden- 
ing to account for the finite MDC width at E^?. Matrix 
elements are included to suppress one side of the disper- 
sion by multiplying the results by a step function with 
a finite width. Ignoring the fact that there is a gap be- 
tween the valence and conduction bands and perform- 
ing an MDC analysis in this region, the fit shows a kink 
in the dispersion (Fig.4(a, b)) and an increase in MDC 
width (Fig.4(c)). The anomalous region where the MDC 
width increases scales with the gap region. This simu- 
lation shows that disregarding the absence of electronic 
states and simply fitting MDCs can sometimes produce 
misleading results. 

It is clear that in order to extract information on both 
the conical dispersion and the presence of a gap, one 
should use a different method. In Fig. 5 we show two 
cuts, one through the K point and one far away from 
the K point. In contrast with Fig 3, the cuts in Fig 5 
are now taken symmetrically with respect to the matrix 
elements. These cuts can be used to measure both sides 
of the cones and extract dispersions without artifacts. 
The use of this geometry is also important because it is 
possible to resolve two peaks in the energy distribution 
curve (EDC) through the K point [23j, a clear signature 
of a gap. 

The dispersions extracted from the EDCs are also 
shown in Fig. 5 (dotted lines in Figs.5(a-b)). In both 
cases we can identify a region of energy where the MDC 
peaks are non-dispersive (Figs.5(c-d)), which coincides 
with the regions of vertical intensity in the upper panels. 
The resemblance between the two cuts is striking and 
once again suggests that a similar phenomenon is the 
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Ed (indicated by broken line) it decreases. This shows 
that near Ejj there are fewer states, a result that is in 
agreement with the gap scenario. When moving away 
from the K point, the angle integrated intensity curve 
(black curve) shows a larger anomalous region, which is 
in agreement with a larger gap. 

In summary, ARPES data taken in single layer epi- 
taxial graphene reveals the opening of a gap at the K 
point. A possible explanation for this gap is that it is 
induced by the breaking of the A-B sublattice symmetry 
[23 which leads to the rehybridization of the valence and 
conduction band states associated with the same Dirac 
point, resulting in a gap at the K point. In the future, 
the possibility of doping graphene to move the Fermi en- 
ergy inside the gap could eventually make graphene a 
semiconductor suitable for electronic applications. 



FIG. 5: (a,b) Dispersions taken along a symmetric direction 
for cuts away and through the K point. (c,d) MDCs in the 
energies as labeled in Fig.5(a,b). (e) Angle integrated inten- 
sity as a function of energy for data in panel a (black curve) 
and panel b (gray curve). 

cause of the vertical region of intensity around the Dirac 
point (marked by the two horizontal arrows). Fig. 5(e) 
shows the integrated intensity for the data of Figs.5(a,b), 
a quantity proportional to the one dimensional density 
of states. Though this one dimensional angle integrated 
intensity does not correspond to the density of states, 
which is integrated over two dimensions, it gives some 
hint whether a gap is present. At high binding energy 
the dispersion shows an almost linear behavior while near 
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